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ABSTRACT 


This report describes the improvements and the modifications of the NASA Aircraft 
Noise Prediction Program (ANOPP) and the Propeller Analysis System (PAS). 

Comparisons of the predictions and the test data are included in the case studies for the flat 
plate model in the Boundary Layer Module, for the effects of applying compressibility 
corrections to the lift and pressure coefficients, for the use of different weight factors in the 
Propeller Performance Module, for the use of the improved retarded time equation solution, 
and for the effect of the number of grids in the Transonic Propeller Noise Module. The 
DNW tunnel test data of a propeller at different angles of attack and the Dowty Rotol data 
are compared with ANOPP predictions. The effect of the number of grids on the Transonic 
Propeller Noise Module predictions and the comparison of ANOPP TPN and DFP-ATP 
codes are studied. In addition to the above impact studies, the transonic propeller noise 
predictions for the SR-7, the UDF front rotor, and the support of the enroute noise test 
program are included. 
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1 . 


Introduction 


The NASA Aircraft Noise Prediction Program (ANOPP) Propeller Analysis System 
(PAS) is a set of computational modules for the prediction of the aerodynamics, 
performance, and noise of subsonic and transonic propellers. The theoretical manual for 
the ANOPP-PAS was published as NASA TM-83199, Part 3 in June 1986. This report 
documents the additions and improvements that have been made to the modules in this 
system. 

Three additional functional modules are available in the ANOPP-PAS for the 
aerodynamic prediction analysis. They are the Improved Blade Shape Module (IBS), the 
Improved Blade Section Aerodynamics Module (IBA), and the Improved Blade Section 
Boundary Layer Module (IBL). These modules are modified versions of functional modules 
RBS, RBA, and BLM in the original system. The IBS module includes a revised blade 
geometry input format and an additional output table for the maximum thickness location 
and trailing edge thickness. The IBA module includes options to use compressibility 
corrections to compute section lift coefficients and local pressure coefficients. The IBL 
module includes an option to use a zero pressure gradient flat plate model to compute the 
boundary layer and a separation of section lift and drag coefficient output tables. Also, the 
boundary-layer thickness is added to the thickness table to match with the Rotor Broadband 
Noise Module (RBN) input table. 

Modifications have also been made to some of the original functional modules in the 
ANOPP-PAS. The Transonic Propeller Noise Module uses a revised method for solving the 
retarded time equation. The Performance Module, the Propeller Loads Module, the 
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Subsonic Propeller Noise Module, and the Transonic Propeller Noise Module have been 
modified to include an option to use the output created by the improved blade aerodynamics 
modules, IBS, IBA, and IBL. A detailed theoretical description of the new functional 
modules and the revisions in each original functional module are presented in Chapter 2. 

Impact studies on the changes to the ANOPP-PAS are documented in Chapter 3. 

The Subsonic Propeller Noise Comparison of the Flat Plate and the Integral Formulation 
Boundary Layer Models presents the thickness and loading noise for each model using test 
case AN5 for the “N” propeller from the FAA DNW Test (reference 1 ). The effects of 
different combinations of compressibility corrections for pressure loading and lift coefficients 
on total noise and the effect of using different weighing factors on the Propeller 
Performance Module are presented using test case AN5 for the “N” propeller from the FAA 
DNW Test. Loading noise comparisons were made using the output of both the original and 
improved blade aerodynamics modules and TPN with and without the new retarded time 
solver. Subsonic propeller noise predictions are made for the R212 propeller (NACA 16 
series blade sections) and compared with the Dowty Rotol data (reference 2). In addition to 
the subsonic propeller noise study, the comparisons of ANOPP predictions and the FAA 
DNW test data are presented for the propellers at different angles of attack. Transonic 
propeller noise predictions are made to study the noise effect of the chosen number of 
grids. The chosen grid combinations are then used for the ANOPP TPN prediction and 
compared with the DFP-ATP (now called ASSPIN) predictions. Further use of the grid 
combinations is made in support of the SR-7 propeller, the UDF front rotor, and the enroute 
noise test program. 

The conclusions from these validation studies are presented in the final chapter. 
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SYMBOLS 


Joukowski transformation parameter, re R 
chord length, re R 
skin friction coefficient 

section lift coefficient 

incompressible lift coefficient 

ambient speed of sound, m/s (ft/s) 

pressure coefficient 
corrected pressure coefficient 

limiting pressure coefficient 

function defining blade surface S 
local Mach number 

ambient Mach number 

blade length measured from axis to tip, m(ft) 

Reynold’s number based on blade length, speed of sound, and kinematic 
viscosity 

leading edge radius, re C 


Reynold’s number 



t time at which noise signal is received by observer, s 

t max maximum thickness, re C 

t te trailing edge bluntness, re C 

free stream velocity, m/s (ft/s) 

x, y surface coordinate, re C 

x t maximum thickness location, re C 

x 2 tr lower trip location, re C 

x 2 tr u upper trip location, re C 

a angle of attack, deg 

8 boundary layer thickness, re C 

8 1 displacement thickness, re C 

5 2 momentum thickness, re C 

y ratio of specific heats 

v kinematic viscosity, m 2 /s (ft 2 /s) 

Q angular velocity of blade, rad/s 

<|> solution of retarded time equation as solved by Newton’s method Q (x - 1) 

\\f blade surface elliptical radial coordinates, rad 

x time at which noise signal is emitted at source position, s 
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spanwise surface coordinate, 0 < ^ < 1 


5i 

chordwise surface coordinate, 0 < 

s leading edge stagnation point, rad 

Subscripts 

i lower surface 

te trailing edge 

tr trip location 

u upper surface 

©o ambient 


5 


2 . 


Functional Modules 


2.1 Improved Blade Shape Module 

INTRODUCTION 

The Improved Blade Shape Module generates a functional representation of the 
airfoil surface using the same methods as the Blade Shape Module (reference 3), with input 
and output table modifications. The improved module incorporates a more concise blade 
section input table and produces an additional output table containing maximum thickness 
location and trailing edge thickness. 


SYMBOLS 

C chord length, re R 

R blade length measured from axis to tip, m (ft) 

t max maximum thickness, re C (r| 2 ) 

t te trailing edge bluntness, re C (rj 2 ) 

x,y surface coordinates in complex z - plane 


x t maximum thickness location, re C (r| 2 ) 

Subscripts: 

t lower surface 

u upper surface 
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INPUT 


This module requires a description of the airfoil surface using the same methods as 
the Blade Shape Module (reference 3). 

Blade Section Table 

The description of the airfoil is entered as a sequence of cross sections defined in a 
Cartesian coordinate system (reference 3). For each cross section, a description array 
containing the spanwise coordinate, leading edge abscissa, leading edge ordinate, chord 
length, leading edge radius, and number of (x,y) pairs on the upper and lower surfaces is 
entered. In the Improved Blade Shape Module, the first four entries in this description array 
are entered normalized with respect to blade length. If the (x,y) coordinates of the upper 
and lower surfaces of a cross section are not identical to the coordinates of the next cross 
section, then the (x,y) coordinates are entered directly after the description array, (i.e. 
upper and lower surface (x,y) coordinates for adjacent identical cross sections are not 
repeated.) 


OUTPUT 

The outputs for this module include all of the outputs produced by the Blade Shape 
Module (reference 3) and an additional table containing the maximum thickness location 
and the trailing edge thickness for each cross section. 
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Blade Thickness Table 


^ array of spanwise coordinates, re R 

x t (^) chordwise location where maximum thickness occurs, re C (^) 
t te (^) thickness of the airfoil cross section at the trailing edge, re C (^) 

METHOD 

Maximum Thickness Location and Trailing Edge Thickness Table 
The maximum thickness of an airfoil section is given by 

^max = (Yu ' ypmax 

The maximum airfoil thickness location x t is the chordwise location where maximum 

thickness occurs. The trailing edge thickness t te is the thickness of the airfoil cross section 
at the trailing edge (i.e. x = 1 ). 


8 



2.2 


Improved Blade Section Aerodynamics Module 


INTRODUCTION 

The Improved Blade Section Aerodynamics Module is a modified version of the 

Blade Section Aerodynamics Module (reference 4). 

The improved module includes the following modifications: 

• The number of Fourier series terms used in evaluating the conformal mapping 
coefficients has been increased. 

• An option to use the Glauert compressibility correction, the Karman-Tsien 
compressibility correction (reference 5), or no correction when computing the pressure 
coefficients has been added. 

• An option to use the Glauert compressibility correction (reference 5) or no 
compressibility correction when computing the lift coefficient has been added. 

• A limiting pressure coefficient check for high negative pressure coefficients based on 
the Carlson and Walkley method has been added (reference 6). 

• Separate tables for the lift coefficients and stagnation points are produced. 

SYMBOLS 

C chord length, re R 

C p pressure coefficient 
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corrected pressure coefficient 



M 

R 


Y 

v 


^2 


limiting pressure coefficient 
incompressible lift coefficient 

section lift coefficient 
incompressible pressure coefficient 
ambient speed of sound, m/s (ft/s) 

ambient Mach number 

local Mach number 
blade length, m (ft) 

Reynold’s number 

Reynold’s number based on blade length, sound speed, and kinematic 
viscosity 

maximum thickness location, re C ( r\ 2 ) 
ratio of specific heats 
kinematic viscosity, m 2 /s (ft 2 /s) 
spanwise surface coordinate, 0 < ^ < 1 
chordwise surface coordinate, 0 < % 2 < 1 
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INPUT 


This module requires the airfoil elliptical coordinates produced by the Improved Blade 
Shape Module (IBS) and arrays of Mach numbers and angles of attack. It requires all of the 
user parameters required by the Blade Section Aerodynamics Module (reference 4) and the 
following additional user parameters. 

C ambient speed of sound, m/s (ft/s) 

oo 

y ratio of specific heats 

R jnf Reynold’s number, based on blade length, speed of sound, and kinematic 

viscosity 


OUTPUT 

This module produces the same output as the Blade Section Aerodynamics Module 
(reference 4) but creates a separate lift coefficient table and separate stagnation point table. 

£ array of spanwise surface coordinates, 0 < ^ < 1 

a array of angles of attack, deg 

M array of Mach numbers 

(^ ,a,M) section lift coefficient 
b >2s (^.a.M) leading edge stagnation point, rad 
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METHOD 


Coefficient of Pressure 

Three options are now available for applying compressibility effects adjustments to 
the pressure coefficients C p found from Bernoulli’s equation (equation 26 in reference 4). 

The first option is to use equation 26 with no compressibility correction. The second option 
applies Glauert’s pressure correction formula 


C 


p - 


Cpo 

Vi - M 2 


( 1 ) 


where M in the equation is set at 0.7 if the actual Mach number is greater than 0.7, 
because the Glauert compressibility correction is not valid at Mach numbers greater than 
0.7. 

The third option applies the Karman-Tsien pressure correction formula (equation 27 in 
reference 4). 


Limiting Pressure Coefficient 

The limiting pressure coefficient is applied to replace extremely high negative 
pressure coefficients on the upper surface of an airfoil. The Reynold’s number based on the 
maximum thickness location is 

R„ = 2 x, C R jn , 

The relation of the limiting pressure coefficient, Mach number and Reynold’s number is 
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IE 


c 


P|im 


-2 

fi 


R x 10’ 6 

n 

R x 10' 6 + 10® 

n 


( 2 ) 


where 

e = 4 - 3M +4(1 - M ) 15 

oo oo 

and 

E = 0.05 + .35 (1 - M J 2 


If the pressure coefficient is negative and less than 2C Pljm then 
C pc = 2C P|im 

The method for calculating the pressure coefficients continues as described for the 
Blade Section Aerodynamics Module (reference 4). 


Coefficient of Lift 

The coefficient of lift for the section is calculated as described in reference 4. There 
are now two options for applying compressibility effects adjustments to the lift coefficients 

C| (equation 25 in reference 4). The first option is to use the lift coefficient with no 

compressibility correction, C| 0 . The second option applies Glauerts compressibility 

correction formula 

r C io 

where M in the equation is set at 0.7 if the actual Mach number is greater than 0.7, 
because the Glauert compressibility correction is not valid at Mach numbers greater than 

0.7. 
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2.3 Improve d Blade Section Boundarv-Laver Module 


INTRODUCTION 

The Improved Blade Section Boundary-Layer Module is a modified version of the Blade 
Section Boundary-Layer Module (reference 7). The improved module computes the two- 
dimensional boundary-layer on airfoil sections using either the integral formulations for the 
boundary-layer thickness (reference 7) or the zero pressure gradient flat plate model 
(reference 8). The user-designated trip locations which locate the transition points are input 
in x coordinates normalized by chord length instead of the elliptical coordinates required by 
the original module. The improved module creates a separate table which contains the 
section drag coefficients. 


SYMBOLS 

b Joukowski transformation parameter, re R 

C chord length, re R 

C f skin friction coefficient 

ambient speed of sound, m/s (ft/s) 

R blade length, m (ft) 

R n Reynold’s number 

r, leading edge radius, re C 

free stream velocity, m/s (ft/s) 
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x chordwise transition point location, re C 

x 2 tr L lower trip location, re C 

x 2 tr u upper trip location, re C 

\j/ blade surface elliptical radial coordinates, rad 

5 boundary-layer thickness, re C 

5 1 displacement thickness, re C 

8 2 momentum thickness, re C 

v kinematic viscosity, m 2 /s (ft 2 /s) 

^ spanwise surface elliptical coordinate, 0 < ^ < 1 

chordwise surface elliptical coordinate, 0 < ^ 1 

Subscripts: 

I lower surface 

tr trip location 

u upper surface 

°o ambient 

INPUT 

This module requires the same input as the Blade Section Boundary-Layer Module 
(reference 7). If the trip locations are required, the x coordinates for the upper surface 
and the lower surface of the trip locations normalized by chord length are input. 
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array of spanwise surface coordinates, 0 < ^ < 1 


Si 

x 2,tr,u (£i) upper trip location, re C 

x 2 tr 7 (£i) lower trip location, re C 


OUTPUT 

This module creates the same output as the Blade Section Boundary- Layer Module 
(reference 7). The section drag coefficients are stored in a separate table. 

^ array of spanwise surface coordinates, 0 < ^ < 1 

a array of angles of attack, deg 

M array of Mach numbers 

C d ,a,M) profile drag coefficient 


METHOD 

Boundary-Layer Equation 

There are two methods that may be used to compute the boundary-layer thickness: 
the integral formulations method described in reference 7 and the zero pressure gradient 
flat plate model. 


Zero Pressure Gradient Flat Plate Model 

For a zero pressure gradient flat plate model, the skin friction coefficients C f increase 
with the distance from the leading edge x (normalized by chord length) and are computed 
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from the Reynold’s number 


C, 


0.0576 



( 4 ) 


or 


C f = 0.0576 (R n .x) 


- 1/5 


At the trailing edge, the turbulent boundary layer thickness 
5(x) = 0.37 1/5 


or 


8<x) = 0.37 (R n x) 


1/5 


the displacement thickness 

8, (x)=5|l 


and the momentum thickness 
S 2 (x) = ^-8(x) 

are computed using the Reynold’s number 

U C 
R = — - 


or 


where 


R „ - R in. MC 


Rinf — " 


C_R 


Trip Locations in x Coordinates Normalized by Chord 
For convenience, the trip locations are input in x coordinates normalized by chord 

and converted to elliptic coordinates % 2 . 

With the input trip location in x coordinates normalized by chord and the assumption 
that the initial elliptical radial coordinates are equal to 0, the elliptic coordinates are 
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computed using 


x = (2b • cosh \j / • cos 2k % 2 + 2b + /2)/C 


( 5 ) 


The new value of the elliptical radial coordinate is then interpolated for the calculated 
elliptic coordinate. The x coordinate is obtained using equation (5) and the whole 
procedure is repeated until the computed x coordinate is equal to the input trip location. 
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2.4 Transonic Propeller Noise Module 


INTRODUCTION 

The Transonic Propeller Noise Module has been modified to include an improved 
retarded time equation solution method developed by Dunn (reference 9). The module also 
includes an option to use the output from the original blade geometry and aerodynamics 
modules, RBS, RBA, and BLM or the improved blade geometry and aerodynamics modules, 
IBS, IBA and IBL. 

SYMBOLS 

function defining blade surface S 
time at which noise signal is received by observer, s 

time at which noise signal is emitted at source position, s 

solution to retarded time equation as solved by Newton’s method, Q. (x-t) 

angular velocity of blade, rad/s 


g 

t 

X 

<l> 

Q 


METHOD 

Roots of the Retarded Time Equation 
The retarded time equation is written as 

g = A<f> 2 + B<t> + C + cos (<]) + D) = 0 (6) 
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where <)> = D(x-t) and only the roots (<j)/Q)<0 are of interest. In order to solve numerically 
for the roots by Newton’s method, the number of roots and a close approximation for each 
root must be determined. 


Step 1 - Define the interval containing all solutions (Bracket the roots): 

Figure 1.1 illustrates that the roots of equation (6) are the points of intersection of 
the parabola A<> 2 + B<{> + C and the sinusoidal curve - cos(<J> + D). Since |cos (<|> + D)| < 1 , 
then the roots of g are bounded by the negative roots of A<|> 2 + B<}> + C = ± 1 , defined to be 
(J) Q and (j)^ . 


Step 2 - Determine the inflection points of g that lie in [<|> 0> (J^J: 

Subintervals [fy, <t> j+1 ] are formed by the endpoints and/or the adjacent inflection 


points in [<|> 0 , <» 1 ]. In each subinterval the curvature of g does not change sign, therefore 
g exhibits parabolic behavior. 


The inflection points ty, <t> i+1 ...<J) N are found by solving 


g" (<» = 2A - cos {<(> + D) = 0, 

I cos-i (2A) - D±2n n 
<b= to < 

I -cos-1 (2A) - D±2m n ’ 0 " 


(n, m = 1, 2, 3 ....) 


( 7 ) 


Step 3 - Determine the number of roots that exist in each subinterval [< 1 ) 1 , <)) j+1 ]: 
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Several tests are made to determine if there are 1 , 2, or no roots in each 


subinterval [fy, <t> i+1 ]. 

(1 ) if g(4>j) * g(<t> i+1 ) = 0 then check the end P oints of the interval MV ^i+il for r00ts 
if g^) = 0, there is a root at 

if g(<|> i+1 ) = 0, then there is a root <f> j+1 

(2) if g(<t>j) • g(<|) i+1 ) <0, then there is exactly one root in [ty, <> i+1 ] 

(3) if g(<})j) • g(<t> i+1 ) >0, then more work has to be done to find roots in 

(a) if g'(<» • g'(<f) i+1 ) > 0, then there are zero roots in [<f>j,<t> i+1 ], if not: 

(b) on [<t>j,<t>j + i]. 9(0) is approximated by a parabola, p(<», with extreme point 

* 

^ext 

(c) find the root of g' in <> ext , by the hybrid Newton/Bisection method 

* 

using ^ext as the initial guess 

(d) if g(<t> ext ) = 0, then there is exactly one double root in [0j,<t>i +1 ] 

(e) if g(<t>j) • g«> GXt ) <0, then there are two roots in [<j> i( <t> i+1 ]: 

(i.e. 1 root in <J> ext ] and 

1 root in [<> ext , <|) i+1 ]) 
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if g(<f>j) * g(0 ext ) > 0 then no roots. 

Step 4 - Determine an initial guess for each root: 

If roots are detected in [<|>j, <t> i+1 ], then g(<j>) is approximated by a parabola p(<j>), in 

the subinterval, and the roots of p(<|>) are used as an initial guess for the roots of g. 

Step 5 - Test each initial guess for stability in Newton’s method: 

Newton’s method can be unstable if ^ is small in some neighborhood of a root. 

d<[) 

If g'(<j>) is less than a given tolerance at the initial estimate of a root then the Bisection 
method is used to further refine the estimate. 

Step 6 - Newton’s method is used to determine the roots: 

The roots of equation (6) are determined through an iterative process, 
a _ a + B<J> n + C + COS (<J) n + D)] 

* n+1 ‘ [2A<f> n + B - sin (<J> n + o)] 

which terminates when |(|) n+1 - <}) n | is smaller than a user-specified convergence criterion or 
if the number of iterations exceeds a limit set by the user. 


22 



2.5 Additional Modifications in ANOPP-PAS Functional Modules 

The Propeller Performance Module (PRP), the Propeller Loading Module (PLD), the 
Subsonic Propeller Noise Module (SPN), and the Propeller Trailing Edge Module (PTE) 
each include an option to use the output from the original blade geometry and 
aerodynamics modules, RBS, RBA, and BLM, or the output from the improved blade 
geometry and aerodynamics modules, IBS, IBA, and IBL. 
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3. 


Impact Studies 


3.1 Impact Study of Boundary Laver Models 

Predictions using the Subsonic Propeller Noise Module for the test case AN5 from 
the FAA DNW test, as described in reference 1 , were made with the following parameters: 


Pitch Angle 

20.8 deg 

RPM 

2700 

Flow Velocity = 

77.0 m/s 

Power = 

184.6 kW 

Thrust 

1907 N 

Attitude Angle = 

0.0 deg 

Flow Temp 

16°C 

Flow Pressure = 

99480 Pa 

Flow Density = 

1.194 kg/m 

Advance Ratio 

0.2680 

Angle of Attack 

1.134 deg 

Power Coefficient = 

0.0490 

Thrust Coefficient 

0.0463 

Helical Tip Mach No.= 

0.8720 


The “N” propeller, used in the FAA DNW test case, is the F8475 D-4 Hartzell, 2-blad- 
ed, 2.03 m diameter propeller with Clark Y airfoil sections, identical to the propeller used on 
the Piper Cherokee Lance aircraft, and has a thin round tip with a thickness ratio of 6.4 per- 
cent at the 3/4 radius position. The two observer coordinates are shown in figure 2.1 follow- 
ing the ANOPP convention and the DNW tunnel convention. 
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The predictions were made using both the flat plate zero pressure gradient boundary- 
layer model and the integral formulation boundary-layer model. The thickness noise and 
loading noise predictions using each model were compared. The predictions at two observ- 
er positions are presented in table 1 . Figures 3.1 and 3.2 show comparisons of flat plate 
and integral formulation data with the DNW data. The results show that the predictions 
using the two models are almost identical and very good in comparison with the DNW data. 
The CPU times required for each of the two models in the Improved Boundary- Layer 
Module are quite different. On MicroVAX, the CPU times required for the execution of the 
Boundary-Layer Module using the integral formulation model is approximately about 30 min- 
utes. Functional module IBL executes 70 times faster when using the flat plate model than 
it executes when using the integral formulation model. The results of this study indicate that 
the zero pressure gradient flat plate model can be used to decrease execution time with lit- 
tle or no loss of accuracy. 
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3.2 Impact Study of the Effects of Applying Compressibility Corrections to 

Pressure and Lift Coefficients 

The Improved Blade Section Aerodynamics Module (IBA) has an option to include 
the effects of Mach number on the lift coefficients using the Glauert theories or to use the lift 
coefficient independent of Mach number. The variation of the pressure coefficients with 
Mach number can be expressed in terms of Glauert or Karman theories, or the pressure 
coefficient can be determined independent of Mach number. For Mach numbers greater 
than 0.7, the Glauert compressibility correction is not valid and the correction for 0.7 Mach 
number is used. 


Predictions for test case AN5 using the “N” propeller from the DNW test were made 
for all possible option combinations. Table 2 provides a summary of the total noise 


predications made for two observer positions using the following option combinations: 


CASE Lift Coefficient 


Pressure Coefficient 


CASE A 
CASE B 
CASEC 
CASE D 
CASE E 
CASE F 


CONSTANT (ICL=0) 
CONSTANT (ICL=0) 
CONSTANT (ICL=0) 

Glauert theory applied (ICL=1) 
Glauert theory applied (ICL=1) 
Glauert theory applied (ICL=1 ) 


CONSTANT (ICP=0) 

Glauert theory applied (ICP=1) 
Karman theory applied (ICP=2) 
CONSTANT (ICP=0) 

Glauert theory applied (ICP=1) 
Karman theory applied (ICP=2) 


The two observer positions are at 60 ° and 90° polar directivity angles as shown in figure 
2.1 . Figure 4.1 and figure 4.2 show the total noise comparisons of the six cases and DNW 
data for the two observers. Figure 4.3 and figure 4.4 present the overall noise of observer 1 
when the Glauert theory is used for the lift coefficients (figure 4.3) and no compressibility 
correction is used for the lift coefficients (figure 4.4). The same for observer 2 is shown in 
figure 4.5 and figure 4.6. 
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The Glauert compressibility correction applied to the lift coefficients provided lower noise 
levels compared to the results from using the lift coefficients independent of the Mach num- 
ber (figures 4.3 and 4.4 for observer 1 ; figures 4.5 and 4.6 for observer 2). For the pressure 
coefficients, the Karman theory provided higher noise compared to the noise results from 
using the Glauert theory. The variation of pressure coefficients with Mach number may be 
expressed by either theory. The Karman correction is more accurate than the Glauert 
correction; however, it is more cumbersome in application. 
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3.3 Impact Stu dy of the Propeller Performance Module 


The Newton iteration method for the induced axial and tangential velocity compo- 
nents in the Propeller Performance Module did not always converge when Glauert com- 
pressibility correction was applied for the lift coefficient. 

The following modifications of the Propeller Performance Module have been accom- 
plished: 

• The interval of the derivatives of the lift and drag coefficients with respect to the 
angle of attack and Mach number is reduced to a smaller interval (from 0.1 to 
0.01 ) for better accuracy. 

• The induced axial velocity component a, and the induced tangential velocity 

component a 2 are initialized for each spanwise station when originally a 1 and 
a 2 were initialized only for the first spanwise station. 

• The weighting factor for the Newton iteration was initially set to 0.5 for the first it- 
eration and set to 1 .0 for the next iteration. Now the weighting factor is set as a 
user parameter and can be varied from 0.1 to 0.8 with a default value of 0.25. 

The effect of using different weighting factors in the functional module is studied 
using the test case AN5 and “N” propeller for observer 1 . Several weighting factors, wt, 
are used for the convergence of the induced velocity. The reference power coefficient, Cp, 

is 0.049. The thrust coefficient C T , power coefficient Cp, efficiency rj , and advanced ratio 

J values are listed in table 3 for comparison. Table 4 presents overall noise in dB for the 
use of different weighting factors. 
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The results show that the effects of the weighting factor are insignificant on the 
convergence of power and thrust coefficients and on the predicted noise levels. Therefore, 
the weighting factor can be varied for the convergence of Newton’s iteration without 
affecting the results. Convergence has been achieved in problems with the Glauert com- 
pressibility correction applied to the lift coefficient. 
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3.4 Impact Study of the Improved Retarded Time Equation Solution 


Predictions were done to study the effects of using the improved retarded time equa- 
tion solution. The pressure coefficients in terms of Glauert and Karman theories were also 
used to provide variable test cases. 

Predictions for the study used the SR-3 propeller which has 4 blades with a 45° tip 
sweep distribution tailored for noise reduction. It also showed the lowest noise level, about 
5 decibels less than the straight-blade SR-2 at Mach 0.8 cruise, and yielded the highest pro- 
pulsive efficiency, 78.7 percent - an improvement of approximately 3 percent over the 
straight-blade SR-2 (reference 10). 

ANOPP coordinates for the two observers located just behind the disk plane in 
meters are as follows: 



0.806 0.00 -0.01 OBSERVER 1 

0.806 0.00 -0.806 OBSERVER 2 

The following prediction test case combinations were used: 

C9?<? Pressure Coefficient Retarded Time Equation Solution 


A Glauert theory applied Old 

B Karman theory applied Old 

C Glauert theory applied New 

D Karman theory applied New 

and the following operating conditions: 


Flight Mach Number, M z = 0.80 
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Helical Tip Mach Number, M h = 0.854 

RPM = 7878.4 

The total noise results for the two observers are shown in table 5. There is less than 
a 0.5 dB difference in the first 6 harmonics comparing case A and case C of observer 1 , a 
1.5 dB difference in the seventh harmonic, and a 4 dB difference in the eighth harmonic. 

For other cases of the two observers, the results from using the old Transonic Propeller 
Noise Module (TPN) and from the modified TPN are identical. The solutions from the im- 
proved retarded time equation solution method have been successfully achieved in all dem- 
onstration problems. 

Newton’s method usually converges in 3 to 5 iterations. The example function is 
shown in figure 5.1 . For the function which is flat as shown in figure 5.2 and has a first de- 
rivative as shown in figure 5.3, more than 10 iterations were required. Therefore, a call to 
the Bisection method was added before a call to Newton's method, and the number of 
permitted iterations was increased to 20 to solve for the flat functions. 
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3.5 A Comparison of ANOPP Predictions with Dowtv Rotol Data 

Predictions were made for the 4-bladed Dowty Rotol R212 propeller having NACA 
16-series blade sections. The Improved Blade Shape Module (IBS), the Improved Blade 
Section Aerodynamic Module (IBA), the Improved Boundary-Layer Module (IBL), the 
Propeller Performance Module (PRP), the Propeller Loading Module (PLD), and the 
Subsonic Propeller Noise Module (SPN) were used for the predictions. 

The propeller is 3.66 meters in diameter and the blade chord, thickness/chord ratio, 


and the blade twist distributions are show in figure 6.1 . Predictions were done for the oper- 
ating conditions as follows: 

Parameter 

Case 1 

Case 2 

Case 3 

RPM 

1250 

1401 

1400 

Flow velocity, m/s 

48 

48 

30 

Power coefficient 

.024 

.036 

.024 

Thrust coefficient 

.023 

.044 

.036 

Helical Tip Mach 

.71 

.794 

.787 

3/4 span blade pitch, deg. 

14.03 

14.03 

8.68 


The near-field microphone is 2.44 m from the propeller hub and the far-field micro- 
phone is 5.49 m from the propeller hub. 

The Dowty Rotol data were from experiments run in the RAE 24-foot anechoic wind 
tunnel (reference 2). Predictions are compared with the experimental data in figure 6.2 for 
case 1 , figure 6.3 for case 2 where the near field and far field microphones are included, 
and figure 6.4 for case 3. Figure 6.2 shows the test case where the propeller helical tip 
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Mach number is moderate and the rpm is low. The first two harmonics show good agree- 
ment, and a difference of 4 dB is shown at higher harmonics. Only the first eight harmonics 
are presented since the data does not extend beyond the frequency of 672 Hz for this case. 
The overall agreement of prediction is good over the range of available experimental data. 
Figure 6.3 shows the case where the helical tip Mach number and rotational speed are high- 
er than case 1 . Agreement is good between experimental data and prediction across all 
harmonics. The third case has a higher propeller speed and a lower power coefficient than 
the first case. Figure 6.4 shows good agreement at lower harmonics and fair agreement at 
higher harmonics. 
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3.6 Impact Study of Propellers at Different Angles of Attack 


Noise predictions were done using the Subsonic Propeller Noise Module (SPN) to 
study the effects of noise for a propeller at different angles of attack. The angles of attack 
are in terms of the propeller shaft inclination and are those from the nacelle axis which were 
tested in the DNW tunnel. 

The “N” propeller was used for this test case. The two microphone positions at 60° 
and 90° polar directivity angles according to ANOPP convention, and the nacelle attitudes 
are shown in Figure 7.1 . 

The operating conditions for the DNW tunnel test cases are as follows: 


Parameter 

9 


Case 3 



RPM 

2700 

2700 

2700 

2700 

2700 

Flow velocity, m/s 

77.3 

76.9 

77.2 

77.0 

77.3 

Power, kW 

152.1 

145.9 

151.6 

149.9 

154.9 

Thrust, N 

1500 

1422 

1432 

1476 

1530 

Nacelle Attitude Angle, deg 

0 

-3.8 

3.6 

-7.4 

7.3 

3/4 span blade pitch, deg 

19.9 

19.9 

19.9 

19.9 

19.9 


The results are compared with the experimental DNW data. For the microphone at 
90° polar directivity angle, figures 7.2 and 7.3 present the total noise at 0°, -3.8°, 3.6°, -7.4°, 
and 7.3° nacelle attitudes. Similarly, figure 8.1 and figure 8.2 present the total noise for the 
60° polar directivity angle microphone at the same above nacelle attitudes. Case 1 is the 
case where the propeller shaft is at 0°. The agreement is good for the first six harmonics 
with the exception of the second harmonic where there is a slight overprediction for both mi- 
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crophones. For the 90° microphone, underprediction appears to be less than 3 dB for the 
higher harmonics. Furthermore, predictions seem to agree better with data for the 90° 
microphone than for the 60° microphone. Cases 2 and 4 are those where the propeller shaft 
is at -3.8° and -7.4°, respectively. In the above two cases for the 90° microphone, with the 
exception of overprediction for the second harmonic, the overall agreement is excellent for 
the first 1 3 harmonics. For higher harmonics, the agreement is still good for the -7.4° 
attitude, and less than 5 dB overpredicted for -3.8° attitude. For both the 60 and 90 
microphones, the prediction and the DNW data are in good agreement. For the cases 
where the propeller shaft is at 3.6° and 7.8°, the DNW data are higher than the predictions 
for all harmonics for both microphone positions. This significant diference may due to the 
non-uniform inflow calculations which is not included in ANOPP . Further study needs to be 
done to determine the possible causes of this underprediction. 
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3 - 7 impact Study of the Effect of the Number of Grids on Transonic Propeller 

Noise Module (TPN1 Predictions 

The pressure time history is affected by the number of propeller blade computation 
grid points and the number of time points per revolution. It is important that a large number 
of propeller grid points near the leading edge be used and enough time points be used to 
give sufficient accuracy of the pressure time history. 

Sixteen combinations of propeller blade computation grids and number of time points 
per revolution were selected for the study as follows: 


Spanwise Stations Chordwise Stations Time Points 


Case 1 

8 

21 

100 

200 

500 

1000 

Case 2 

12 

31 

100 

200 

500 

1000 

Case 3 

19 

41 

100 

200 

500 

1000 

Case 4 

25 

51 

100 

200 

500 

1000 


The two observer coordinates were chosen as follows: 



x i 

X 2 

X 3 

Observer 1 

3.06 

0.0 

0.833 

Observer 2 

3.06 

0.0 

0.000 


The SR-7 propeller, one of the propellers used in the Advanced Turboprop Project 
(ATP), was used for these preliminary studies. This propeller is similar to the SR-3 but has 
less sweep (41 ° versus 45°). The flight conditions for this study were a Mach number of 
0.8, a flight altitude of 35000 feet, and a propeller speed of 1692 RPM. 
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In order to determine the minimum required number of time points for accuracy, the 
first three harmonics of each observer with the four grid cases for each time point number 
were compared. Figures 9.1 to 9.3 and figures 10.1 to 10.3 show that for the time points of 
500 and 1000, the difference in the total noise is insignificant. For the four cases, especially 
cases 1 and 2 where there are less propeller grid points, the same noise levels are pro- 
duced. Therefore, the minimum required number of time points for the analysis was chosen 

to be 500. 

A large number of propeller grid points in the regions where the blade is highly 
curved is recommended. In the combination grid cases, the difference in amplitude 
between cases 3 and 4 is fairly small compared to that between cases 1 and 2. Therefore, 
case 3, where there are 19 spanwise and 41 chordwise grid points, was selected as having 
the optimum number of grid points to save computation time. 
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3.8 The Comparison of ANOPP TPN and ASSPIN Codes 


The ASSPIN code is the latest theoretical acoustic result for the transonic/supersonic 
sources developed by Mark H. Dunn, F. Farassat, and Sharon L. Padula (reference 9). For 
the noise comparison of the two codes, the same aerodynamic loads generated from Adam- 
czyk’s Euler code are used for input to ANOP TPN code. The propeller grid combination of 
19 grid points for the span and 41 grid points for the chord with 500 time points was used for 
ANOPP TPN. 

A comparison of the ANOPP TPN and the ASSPIN codes was done for the SR-7 at 
1=1 and k=1-20 grid points (figure 11.1). The flight conditions for this study are a Mach 
number of 0.8, a flight altitude of 35000 feet, and a propeller speed of 1692 RPM. The first 
harmonic of the ANOPP TPN and ASSPIN (figure 1 2.1 ) shows that the maximum difference 
of 1 .5 dB is in the area behind the propeller plane. However, in figure 1 2.2, for the second 
harmonic, the noise from ASSPIN prediction is about 2 dB higher than the noise from 
ANOPP TPN in the area in front the propeller plane. The opposite effect is shown in the 
area behind the propeller. For the third harmonic, ASSPIN produces almost 5dB higher 
noise in the aft area of the propeller (figure 1 2.3). The noise difference seems to be greater 
at higher harmonics. Since there is not a great difference in the sound pressure level be- 
tween ANOPP TPN and ASSPIN especially for the first harmonic, either ANOPP TPN or 
ASSPIN can be chosen for transonic noise prediction. Further study should be done for the 
noise comparison at higher harmonics. 
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3.9 Impact Study of the Improved TPN Using the SR-7 on the ATP and the UDF 
Front Rotor on the 727 

The pressure amplitude and phase data for the exterior noise levels were computed 
using the ANOPP predictions. The area of interest was k=5 to 20 and 1=1 to 10 (figure 
1 1.1) for the SR-7. The predictions were used for input to the Propeller Aircraft Interior 
Noise (PAIN) program. A propeller grid having 19 grid points spanwise and 41 grid points 
chordwise combined with 500 time points was used for computation. 

The grid is specified by coordinates x v x 2 , and x 3 as shown in figure 11.1. The co- 
ordinate positions are then given by the following relations for k=1 ,20 and 1=1 ,10: 
x 1 = (r p - a) = a {1 - cos [(I -1 ) tc/1 8]} 

x 2 = -a sin [(1-1) tc/ 18] 
x 3 = (rca/l 8) [k p - k] 
where 

k p is the k index of the grid point (k, I) = (k p , 1), that is, where the coordinate x 1 
penetrates the fuselage 

r p is the radial distance from the center of the fuselage cylinder to the axis of 
rotation of the propeller 

a is the radius of the fuselage 
k p = 8 

r p - 4.2561 5 m 
a = 1.1938 m 

The operating conditions were: 


39 




Flight 

Mach Number 

Flight 

Altitude(ft) 

Propeller 

RPM 

Power 

Coefficient 

0.8 

35,000 

1692 

1.787 

0.8 

25,000 

1777 

1.560 

0.5 

15,000 

1692 

1.33 


The refined grid, which is the mid-point of the space between the two grid points, is 
again used in the prediction for the first flight case. 

The same grid point combination is used for the prediction of the UDF front rotor with 
the following flight conditions: 


Flight 

Flight 

Propeller 

Thrust 

Mach Number 

Altitudefft) 

RPM 

M 

0.8 

35,000 

1267.5 

2100 

0.72 

35,000 

1275 

2100 


The coordinate positions (figure 14.1) are given by the following relations for k=1 to 
16 and 1=1 to 10: 

x 1 = r p - a cos[(l - 1 ) Jt/1 8] 1 - cot a 1 + (k p ' ccrt a 
x 2 =a sin[(l - 1 ) jc/1 8] 1 - -£■ cot a 1 + ^ (k p - k| cot a 
x 3 =( 7 ta/ 1 8) [k p - k 1 - cot a 

a = cos' 1 (a/s 2 ) 

where 

r p is the radial distance from the center of the fuselage cylinder to the axis 
of rotation of the propulser 

a is the radius of the fuselage 
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s 2 is the meridian distance to the apex of the frustum 

r p = 3.3782 m 
a= 1.880 m 
z p = 3.91 m 

s 2 = 9.40 m 

The first harmonics for the 3 cases of the SR-7 at L = 1 are shown in figure 1 3.1 . 
Case 1 , with a power coefficient greater than that of case 2, provides higher amplitude. 

Case 3, which is the subsonic case, has the lowest amplitude compared to cases 1 and 2. 
The overall predictions for case 1 (figure 13.2), case 2 (figure 13.3), and case 3 (figure 13.4) 
show that moving from L=1 (closest to the propeller) to L = 10 (further away from the 
propeller), result in amplitude decreases in front the propeller plane. However, in the area 
aft of the propeller (starting at k=12), the amplitude increases going from L = 1 to 10. 

For the UDF (figure 15.1) at longitudinal positions k= 1 to 9, the amplitude is low 
at L = 1 and progressively increases at L = 10. The trend is opposite from position k = 9 
to the plane of the propeller and continues aft to k=13. Aft of k = 13 the trend with 
changing L returns to that found forward of k = 9. Case 2 (figure 15.2) has the same 
trends with L as seen in figure 15.1 but shifted to different longitudinal positions. 

The amplitude and phase trends seen in figures 15.1 and 15.2 have an impact on 
results from the Interior Noise Prediction Program (PAIN). 
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3.10 Impact Study from the Support of the Enroute Noise Test Program 


Predictions of the PTA propfan noise to be used as input to the ray tracing propaga- 
tion code were made. These predictions used the same propeller computation grid and the 
same number of time points per revolution as were used for the predictions of the PTA 
noise for the Propeller Interior Noise Program (PAIN) study. 

The propeller speed is 1697.65 RPM for all of the following flight conditions: 


Case 

Altitude (ft) 
(Above Sea Level) 

Mach Number 

Power Coefficient 

1 

7,000 

0.5 

1.28 

2 

14,000 

0.5 

1.34 

3 

20,000 

0.5 

1.37 

4 

20,000 

0.7 

1.55 

5 

35,000 

0.7 

1.48 


An atmospheric absorption table was built (using the ATM and ABS modules) with a 
reference altitude of 5,000 feet and 70% relative humidity for the White Sands standard day 
meteorological conditions. The flight altitudes were therefore set relative to the White 
Sands test level. 

For each flight condition, the PTA propfan predictions were done for 25 observer po- 
sitions starting at a polar directivity angle of 30° with 5° intervals and at a source-to-observer 
distance of 500 ft (figure 16.1). 

The first and the second harmonics of the eight requested harmonics were examined 
for the above five cases as shown in figures 17.1 to 17.5. Cases 1, 2, and 3, where the 
flight Mach number was 0.5, provide smooth curves compared to cases 4 and 5, where the 
flight Mach number was 0.7. First harmonic noise comparisons for each of the five flight 
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conditions are presented in figure 18.1 and those for the second harmonic are presented in 
figure 18.2. For the first harmonic, the first three cases show very similar noise levels 
(figures 17.1, 17.2, and 17.3). 
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4. Example Problem 

Usually, the predictions are made in two steps for the use of the subsonic propeller or 
in three steps for the use of the transonic propeller. If the computation requires a large 
number of chordwise and spanwise stations for the airfoil sections, a large number of grid 
points or observers, and a large number of time points per revolution; then the prediction 
should be made in three steps. 

The first step is to build the blade library from the execution of IBS, IBA, and IBL as 
illustrated below: 


CREATE GRID $ 

UPDATE NEWU=GRID SOURCE=* $ 

-ADDR OLDM=* NEWM=PSI FORMAT=4H*RS$ MNR=1 $ 

0 . 0 $ 

-ADDR OLDM=* NEWM=XI2 FORMAT=4H*RS$ MNR=1 $ 

0.0 0.0417 0.0833 0.1250 0.1667 0.2083 0.2500 
0.2917 0.3333 0.37509 0.41667 0.45804 0.50000 
0.54167 0.58333 0.62500 0.66667 0.70833 0.75000 
0.79167 0.83333 0.87500 0.91667 0.95833 1.00000 $ 
END* $ 

PARAM IUNIT=2HSI $ 

CREATE GEOM $ 

UPDATE NEWU=GEOM SOURCE=* $ 

-ADDR NEWM=BLADE OLDM=* FORMAT=0 $ 

10 $ 

.3 0324 * . 1762 ' .340 6.85E-3 79.6 24 23 $ 

1.000000 0.000000 $ 

.990000 .002020 $ 

.980000 .004044 $ 


.4 


.980000 .000660 $ 

.990000 .000356 $ 

1.000000 0.000000 $ 

-.0571 .2159 .374 2.50E-3 75.2 24 23 $ 

1.000000 0.000000 $ 

.990000 .001677 $ 

.980000 .003347 $ 
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.980000 .001357 $ 

.990000 .000699 $ 

1.000000 0.000000 $ 

.5 -.0806 .2350 .400 9.42E-4 71.0 50 49 $ 

1.000000 .000000 $ 

.990000 .002505 $ 

.980000 .004825 $ 


.980000 -.000152 $ 

.990000 -.000144 $ 

1.000000 -.000000 $ 

.6 -.0874 .2054 .4 5.63E-4 66.9 50 49 $ 

1.000000 .000000 $ 

.990000 .002210 $ 

.980000 .004265 $ 


.980000 .000407 $ 

.990000 .000150 $ 

1.000000 -.000000 $ 

.7 -.0730 .1438 .376 3.29E-4 63.1 50 49$ 

1.000000 .000000 $ 

.990000 .001972 $ 

.980000 .003814 $ 


.980000 .000855 $ 

.990000 .000386 $ 

1.000000 -.000000 $ 

.8 -.0350 .0588 .332 2.36E-4 59.2 50 49 $ 

1.000000 .000000 $ 

.990000 .001852 $ 

.980000 .003588 $ 


.980000 .001080 $ 


45 



.990000 .000506 $ 

1.000000 -.000000 $ 

.85 -.0066 .0106 .304 2.15E-4 57.9 50 49 $ 

.9 .0312 -.0468 .258 2.15E-4 56.3 50 49 $ 

1.000000 .000000 $ 

.990000 .001822 $ 

.980000 .003531 $ 


.980000 .001136 $ 

.990000 .000536 $ 

1.000000 -.000000 $ 

.95 .0748 -.1063 .204 1.95E-4 54.9 50 49 $ 

.99 .1143 -.1568 .154 1.95E-4 53.9 50 49 $ 

1.000000 .000000 $ 

.990000 .001792 $ 

.980000 .003474 $ 


.980000 .001192 $ 

.990000 .000566 $ 

1.000000 -.000000 $ 

END* $ 

EXECUTE IBS $ 

PARAM IPRINT=3 NORDER=4 $ 

PARAM ICL=0 $ NO COMPRESSIBILITY CORRECTION FOR 

$ LIFT COEFFICIENT 

PARAM ICP=2 $ KARMAN COMPRESSIBILITY CORRECTION 

$ FOR PRESSURE COEFFICIENT 
CREATE I BA $ 

UPDATE NEWU=IBA SOURCE=* $ 

-ADDR NEWM=MACH OLDM=* FORMAT=4H*RS$ MNR=1$ 

.1 .3 .5 .7 .9 $ 

-ADDR NEWM=ALPHA OLDM=* FORMAT=4H*RS$ MNR=1 $ 

-1. 1. 3. $ 

END* $ 

PARAM CA = 309.7073 $ SPEED OF SOUND IN M/S 
PARAM VNU = .000028026 $ VISCOSITY IN M**2/S 

PARAM B = 1.3716 $ BLADE LENGTH IN M 

EVALUATE RINF = CA*B/VNU $ 

EXECUTE IBA $ 

$ 

$ 

EXECUTE IBL $ 

UNLOAD /SR3LIB1/ $ 

PARAM ICL=1 ICP=2 $ GLAUERT CORRECTION FOR LIFT, 

$ KARMAN CORRECTION FOR PRESSURE 
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EXECUTE I BA $ 

EXECUTE IBL $ 

UNLOAD / SR3LIB2 / $ 
ENDCS $ 


• The second step is to use the libraries built in the first step to build the Loads 


library for the noise calculation. 


LOAD/SR3LIB1 / $ LOAD IN THE GEOM LIBRARY BUILT IN THE FIRST STEP 
PARAM EPSILON=0 . 0 1 


$ 

$ 

$ 

PARAM OPTION=l, UNIFORM= . TRUE . 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

$ 

PARAM IMPROV = .TRUE. 


$ ERROR CRITERION FOR 
STOPPING THE ITERATION 
ON THE INDUCED ANGULAR 
AND AXIAL VELOCITIES 
$ OPTION, METHODOLOGY 
OPTION 

= 0: BLADE ELEMENT- 

MOMENTUM THEORY 
= 1 : PRANDTL TIP RELIEF 

CORRECTION 

UNIFORM, INFLOW OPTION 
= TRUE : UNIFORM INFLOW 
SPECIFIED BY MZ 
= FALSE: RAD I ALLY VARYING 
INFLOW SPECIFIED 
BY RBF (FLOW) 


$ FLAG TO USE THE IMPROVED VERSION 
$ OF ANOPP 


CREATE GRIDA $ 

UPDATE NEWU=GRIDA SOURCE=* $ 

-ADDR NEWM=PSI OLDM=* FORMAT=4H*RS$ 

0.0 $ 


)R NEWM= 

XII OLDM=* FORMAT=4H*RS$ 

MNR=1 $ 



.300 

.3380 

.4334 

.4669 




.5004 . 

5980 . 

6956 . 

7337 .7718 

.8060 



.8402 . 

8706 . 

9010 . 

9277 .9544 

.9772 

99 $ 


)R NEWM= 

■XI 2 OLDM=* FORMAT=4H*RS$ MNR=1 $ 



.0000 

.0933 

.1602 

.2337 .2836 

.3225 

.3536 

.3789 

.3940 

.4002 

.4186 

.4267 .4348 

.4421 

.4494 

.4562 

.4629 

.4693 

.4757 

.4819 .4880 

.4940 

.5000 

.5060 

.5120 

.5182 

.5243 

.5307 .5371 

.5439 

.5506 

.5652 

.5814 

. 6211 

. 6464 

.6775 .7164 

.7663 

.8302 

.9096 

1.000 $ 







END* $ 







PARAM B 

=1.3716, IUNITS=2HSI IPRINL 

=3 $ 




PARAM ALPHAP = 0 


PARAM 

PARAM 


IDPDT 

BETA75 


= 0 
= 60 
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$ SET PROPELLER ANGLE OF 
$ ATTACK IN DEGREES 
$ PROPELLER LOADING IS STEADY 
$ INITIAL GUESS FOR PROPELLER 
$ 3/4 SPAN 
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$ PITCH ANGLE IN DEGREES 
PARAM RPM =1777. $ PROPELLER RPM 

PARAM NBLADE = 4 $ NUMBER OF PROPELLER BLADES 

PARAM ORIG =61.15 $ BLADE TWIST FROM ROOT TO 3/4 

$ SPAN 

EVALUATE RHOA= . 14 961 6E-2 *32 . 1726*1 6 . 02 $ DENSITY IN KG/M**3 
PARAM I UNITS = 2HSI $ SI UNITS 

PARAM IPRINT =3 $ REQUEST INPUT AND OUTPUT 

$ PRINT 

$ 

$ EVALUATE CONTROL STATEMENTS ARE USED TO COMPUTE 
$ ADDITIONAL REQUIRED QUANTITIES 
$ 

$ 


EVALUATE 

RPS 

= RPM / 

60. 



$ 

COMPUTE REVOLUTIONS PER SECOND 

PARAM 

PI 

= 3.1415926 



$ 

SET VALUE OF PI 

PARAM CPREF = 1 

.56 $ 


EVALUATE 

BETA 

= BETA75 

- ORIG 



$ 

COMPUTE ROOT PITCH IN DEGREES 

EVALUATE 

THETAR 

= BETA * 

PI / 180. 



$ 

CONVERT ROOT PITCH TO RADIANS 

EVALUATE 

ALPHAP 

= ALPHAP 

* PI / 180. 



$ 

CONVERT PROPELLER ANGLE OF 



$ 

ATTACK TO RADIANS 

PARAM 

CA = 309.70728 




$ 

SPEED OF SOUND IN M/S 

PARAM 

VF 

= 247.766 



$ 

VELOCITY IN M/SEC 

EVALUATE 

MZ 

= VF / CA 



$ 

COMPUTE FORWARD MACH NUMBER 

EVALUATE 

OMEGA 

=2. * PI * RPS 



$ 

COMPUTE ANGULAR VELOCITY 

EVALUATE 

MACHRF 

= B * OMEGA / CA 



$ 

COMPUTE ROTATIONAL TIP MACH 



$ 

NUMBER 


$ 

$ ADDITIONAL OUTPUT CONTROL PARAMETERS ARE 

$ 

$ 

$ 

$ THE 3/4 SPAN BLADE PITCH MUST BE ADJUSTED SO THAT THE 
$ COMPUTED POWER COEFFICIENT MATCHES THE MEASURED 
$ POWER. THIS REQUIRE AN ITERATIVE SOLUTION OF THE 
$ PROPELLER PERFORMANCE (PRP) MODULE. THE SECANT 
$ METHOD IS USED TO FIND THE ROOT TO THE EQUATION F(Z) = 
$ CPREF - CP. CONVERGENCE IS ASSUMED COMPLETE WHEN THE 
$ COMPUTED VALUE IS WITHIN ONE PERCENT OF THE MEASURE 
$ VALUE. 

$ 

$ 
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PARAM Z1 = THETAR $ 

EXECUTE PRP GRID=GRIDA $ 

EVALUATE FZl = CPREF - CP $ 

EVALUATE THETAR = THETAR + PI / 180. $ 

PARAM Z2 = THETAR $ 

PARAM COUNT = 1 $ 

LABI CONTINUE $ 

EXECUTE PRP GRID=GRIDA $ 

EVALUATE FZ2 = CPREF - CP $ 

EVALUATE DIFF = FZ2 / CPREF $ 

EVALUATE DIFF=ABS (DIFF) $ 

EVALUATE COUNT = COUNT + 1 $ 

IF ( DIFF .LT. 0.01 ) GOTO LAB 2 $ 

EVALUATE Z = Z2 - FZ2 * ( Z2 - Z1 ) / ( FZ2 - FZl ) $ 

PARAM Z1 = Z2 $ 

PARAM Z2 = Z $ 

PARAM FZl = FZ2 $ 

PARAM THETAR = Z $ 

EVALUATE COUNT = COUNT + 1 $ 

IF ( COUNT .GT. 10 ) GOTO LAB 3 $ 

GOTO LABI $ 

LAB 2 CONTINUE $ 

EXECUTE PLD GRID=GRIDA $ 

LAB 3 CONTINUE $ 

UPLIST $ 

UNLOAD /PRPLIB1/ PLD $ 


• Finally, the third step is where the noise is calculated. The blade geometry library 
which is built in the first step and the load library which is built in the second step are used 
in this step as shown: 


LOAD / SR3LIB1 / $ 

LOAD /PRPLIB1 / $ 

PARAM IMPROV= . TRUE . $ 

EVALUATE Z = 25000 . *0 . 3048 $ 

EVALUATE RHOA = .106626E-2 * 32.1726 * 16.02 $ 

PARAM PI=3. 1415926, REV=1777 . , IPRINT=3 $ 

EVALUATE CA=1016.10 * .3048 $ 

EVALUATE OMEGA=2 . *PI * REV / 60 . $ 

PARAM MZ = 0.8 $ 

EVALUATE B=4 . 5 * 0.3048 $ 

EVALUATE MACHRF = B * OMEGA / CA $ 

PARAM THETAR=0 . 0 , EPSILON=0 . 00 1 , NBLADE=4 $ 

CREATE TPN $ 

UPDATE NEWU=OBSERV SOURCE=* $ 

-ADDR NEWM=COORD OLDM=* FORMAT=4H3RS$ $ 
3.34 -.77 1.46 $ 

END* $ 


$ 
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{/>{/></> 


BLADE SECTIONS USE A GRID WITH EXTRA POINTS NEAR 
THE LEADING EDGE (IE WHERE XI2 = PI ) 


$ 


$ 


CREATE GRIDA$ 

UPDATE NEWU=GRIDA SOURCE^* $ 

-ADDR OLDM=GRID (XII) $ 

-ADDR OLDM=GRID (PSI) $ 

-ADDR NEWM=XI2 OLDM=* FORMAT=4H*RS$ 

0.0000 0.0933 0.1602 0.2337 0.2836 0.3225 

0.4002 0.4186 0.4348 0.4494 0.4629 0.4757 

0.5120 0.5243 0.5371 0.5506 0.5652 0.5814 

0.6464 0.6775 0.7164 0.7663 0.8302 0.9096 

END* $ 


MNR=1 $ 

0.3536 0.3789 
0.4880 0.5000 
0.5998 0.6211 
1.0000 $ 


PARAM PSI0=0., NT=200 $ 

PARAM NTIME = 500 $ 

PARAM NHARM=8 , IPRES=1 $ 
PARAM TRANS=0 . 98, FRACDT=1 . 0 $ 


EXECUTE TPN GRID=GRIDA $ 


If the job is small and there is no need to reuse the blade geometry library and load 
library then the second and the third step can be combined, or all three steps can be 
combined into one job. 
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5. Conclusions 

The modified version of ANOPP contains three additional modules which are 
improved versions of the Blade Shape Module (IBS), the Rotating Blade Aerodynamic 
Module (IBA), and the Boundary-Layer Module (IBL). In addition, the Propeller 
Performance Module, the Propeller Loading Module, and the Subsonic Propeller Noise 

Module have been modified and improved. 

For the Blade Shape Module, there is no repetition for airfoil sections with the same 

coordinates. For the Rotating Blade Aerodynamic Module, more options are available for 
the compressibility correction of the lift and the pressure coefficients. For the Boundary- 
Layer Module, the addition of the flat plate model option provides no loss in accuracy and 
less amount of computation time. Due to the modification of the Propeller Performance 
Module, the convergence of the Newton method in the Propeller Performance Module is 
fairly rapid, including the case where the Glauert theory is applied for the lift coefficient. As 
a result, the Propeller Loading Module, the Subsonic Propeller Noise Module, and the 
Transonic Propeller Noise Module with the changed retarded time solver work properly now. 

Besides having more options, the improved version of ANOPP can also give the 
same results as the original version if the compressibility correction of the lift coefficient is 

set to 0, and that of the pressure coefficient is set to 2. 

The study comparing the predictions of the improved version of ANOPP, the old 
version of ANOPP, and the available test data shows good agreement between predictions 
and test data. Therefore, the improved version of ANOPP should be used to save 
execution time when necessary with no loss in accuracy. 
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Table 1. Subsonic Propeller Noise Comparison of the Zero Pressure Gradient Flat Plate and the 
Integral Formulations for the Boundary-Layer of Thickness Models 

OBSERVER 1 - 0 = 90°, <}> = 0°, r = 4 m 
OBSERVER 2- 9 = 60°, <f> = 0°, r = 4m 



OBSERVER NUMBER 1 

OBSERVER NUMBER 2 


FLAT PLATE 

NON-FLAT PLATE 

FLAT PLATE 

NON-FLAT PLATE 

HARMONIC 

THICKNESS 

LOADING 

THICKNESS 

LOADING 

THICKNESS 

LOADING 

THICKNESS 

LOADING 

NUMBER 

NOISE dB 

NOISE dB 

NOISE dB 

NOISE dB 

NOISE dB 

NOISE dB 

NOISE dB 

NOISE dB 

1 

104.02 

115.23 

104.02 

114.98 

101.10 

110.68 

101.10 

110.89 

2 

108.40 

114.51 

108.40 

114.25 

103.41 

109.86 

103.41 

110.06 

3 

109.79 

112.96 

109.79 

112.68 

102.82 

107.96 

102.82 

108.13 

4 

110.19 

111.24 

110.19 

110.96 

101.32 

105.53 

101.32 

105.69 

5 

110.14 

109.49 

110.14 

109.21 

99.40 

102.82 

99.40 

102.96 

6 

109.81 

107.75 

109.81 

107.47 

97.23 

99.93 

97.23 

100.06 

7 

109.28 

106.02 

109.28 

105.75 

94.87 

96.92 

94.87 

97.05 

8 

108.60 

104.30 

108.60 

104.03 

92.38 

93.83 

92.38 

93.94 

9 

107.80 

102.59 

107.80 

102.33 

89.77 

90.67 

89.77 

90.77 

10 

106.90 

100.88 

106.90 

100.64 

87.06 

87.45 

87.06 

87.55 

1 1 

105.91 

99.17 

105.91 

98.95 

84.26 

84.19 

84.26 

84.27 

12 

104.84 

97.47 

104.84 

97.26 

81.37 

80.88 

81.37 

80.95 

13 

103.69 

95.76 

103.69 

95.58 

78.40 

77.53 

78.40 

77.59 

14 

102.46 

94.06 

102.46 

93.90 

75.35 

74.15 

75.35 

74.19 

15 

101.17 

92.35 

101.17 

92.22 

72.22 

70.73 

72.22 

70.77 

16 

99.81 

90.64 

99.81 

90.54 

69.01 

67.29 

69.01 

67.31 

17 

98.39 

88.93 

98.39 

88.87 

65.72 

63.82 

65.72 

63.84 

18 

96.89 

87.22 

96.89 

87.19 

62.35 

60.34 

62.35 

60.34 

19 

95.33 

85.50 

95.33 

85.51 

58.89 

56.85 

58.89 

56.84 

20 

93.69 

83.79 

93.69 

83.82 

55.35 

53.35 

55.35 

53.34 
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Table 2. Subsonic Propeller Noise Comparison from the use different theories 
for the computations of Lift and Pressure Coefficients: Karman-Tsien (ICL=2, 
ICP=2), Glauert (ICL=1 ,ICP=1), and Independent of Mach number (ICL=0, 


ICP 



ICL = 0 
ICP = 0 
CASE A 

ICL = 0 
ICP = 0 
CASE B 

ICL = 0 
ICP = 2 
CASEC 

HARMONIC 

NUMBER 

OBSERVER 1 9 = 90° 

, <ho°, 

1 

112.150 

114.210 

115.430 

2 

112.370 

114.070 

115.390 

3 

112.020 

113.330 

114.620 

4 

111.570 

112.540 

113.710 

5 

111.030 

111.750 

112.770 

6 

110.410 

110.940 

111.810 

7 

109.690 

110.100 

110.830 

8 

108.890 

109.200 

109.800 

9 

108.000 

108.240 

108.730 

10 

107.040 

107.220 

107.610 

11 

105.990 

106.130 

106.440 

12 

104.880 

104.980 

105.210 

13 

103.690 

103.770 

103.930 

14 

102.440 

102.490 

102.600 

15 

101.120 

101.150 

101.200 

16 

99.740 

99.750 

99.750 

17 

93.290 

98.280 

98.230 

18 

96.770 

96.750 

96.650 

19 

95.180 

95.140 

95.000 

20 

93.510 

93.470 

93.270 


0 ). 


ICL = 1 
ICP = 0 
CASE D 

ICL = 1 
ICP=1 
CASEE 

ICL = 1 
ICP = 2 
CASE F 

DNW 

Data 

4m 




110.400 

112.330 

113.380 

115.0 

111.170 

112.610 

113.670 

113.1 

111.230 

112.240 

113.190 

115.1 

111.030 

111.730 

112.530 

113.5 

110.660 

111.150 

111.800 

113.2 

110.150 

110.500 

111.010 

111.2 

109.510 

109.770 

110.150 

113.0 

108.760 

108.950 

109.230 

110.0 

107.900 

108.050 

108.240 

109.9 

106.960 

107.070 

107.190 

109.4 

105.940 

106.010 

106.070 

107.3 

104.830 

104.890 

104.890 

106.3 

103.660 

103.690 

103.650 

106.2 

102.410 

102.430 

102.340 

102.7 

101.100 

101.100 

100.960 

103.1 

99.710 

99.700 

99.530 

102.0 

98.260 

98.240 

98.020 

99.5 

96.740 

96.700 

96.440 

97.8 

95.140 

95.090 

94.790 

97.3 

93.480 

93.410 

93.060 

95.7 



CASE A 

CASE B 

CASEC 


OBSERVER 2 0 = 

60°, <|> = 0 o , 

1 

107.350 

109.170 

111.850 

2 

107.210 

108.690 

111.250 

3 

105.650 

106.940 

109.380 

4 

103.530 

104.660 

107.000 

5 

101.160 

102.140 

104.340 

6 

98.620 

99.480 

101.520 

7 

95.970 

96.710 

98.590 

8 

93.230 

93.860 

95.570 

9 

90.410 

90.940 

92.480 

10 

87.510 

87.960 

89.320 

11 

84.550 

84.910 

86.100 

12 

81.520 

81.810 

82.810 

13 

78.410 

78.640 

79.460 

14 

75.240 

75.400 

76.050 

15 

72.000 

72.100 

72.570 

16 

68.680 

68.720 

69.010 

17 

65.290 

65.270 

65.380 

18 

61.810 

61 .730 

61 .660 

19 

58.240 

58.110 

57.850 

20 

54.580 

54.390 

53.930 


DNW 


CASE D 

CASE E 

CASE F 

Data 

= 4m 




106.180 

107.730 

109.910 

111.9 

106.410 

107.600 

109.550 

108.2 

105.100 

106.100 

107.890 

109.0 

103.180 

104.050 

105.710 

106.2 

100.980 

101.740 

103.260 

104.4 

98.590 

99.270 

100.650 

100.0 

96.070 

96.670 

97.930 

101.1 

93.440 

93.980 

95.110 

100.0 

90.710 

91.200 

92.230 

93.9 

87.910 

88.350 

89.270 

94.4 

85.020 

85.420 

86.240 

91.9 

82.060 

82.430 

83.160 

85.5 

79.030 

79.360 

80.010 

85.4 

75.930 

76.230 

76.810 

82.0 

72.750 

73.030 

73.530 

78.1 

69.500 

69.760 

70.200 

76.2 

66.180 

66.410 

66.790 

80.5 

62.770 

62.990 

63.310 

72.6 

59.280 

59.480 

59.760 

65.4 

55.710 

55.890 

56.130 

70.3 
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Table 3. Power Coefficients C pp Thrust Coefficients C T , Advanced 
Ratio J, and Efficiency t| at Variable Weighting Factors 


wt. 

0.25 

0.3 

0.4 

mm 

0.6 

0.7 

0.8 

Co 

0.0490039 

0.0490023 

0.0490006 

0.0490060 

0.0490076 

0.0490083 

0.0490083 

Ct 

0.0496741 

0.0496339 

0.0495821 

0.0495703 

0.0495497 

0.0495405 

0.0495358 

j 

0.843905 

0.843905 

0.843905 

0.843905 

0.843905 

0.843905 

0.843905 

n 

0.855446 

0.854782 

0.853920 

0.853623 

0.853240 

0.853069 

0.852989 

No. of 
Iterations 

3 

3 

3 

3 

3 

3 

3 


Table 4. Overall Noise Comparison for the Use of Variable Weighting Factors 


Wt = .25 

Wt = .3 

Wt = .4 

Wt.= .5 


Wt = .7 

Wt= ,8 

115.38 

115.39 

115.40 

115.40 

115.41 

115.41 

115.41 

115.36 

115.36 

115.37 

115.37 

115.38 

115.38 

115.38 

114.59 

11 4^59 

114.60 

114.60 

114.60 

114.60 

114.60 

113.68 

113.68 

113.69 

113.69 

113.70 

113.70 

113.70 

112.75 

112.75 

112.76 

112.76 

112.76 

112.76 

112.76 

111.79 

111.79 

111.80 

111.80 

111.80 

111.81 

111.81 

110.81 

110.81 

110.81 

110.81 

110.82 

110.82 

110.82 

109.78 

109.78 

109.79 

109.79 

109.79 

109.79 

109.79 
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Table 5. Transonic Propeller Noise (dB) Comparison of the Old and Improved 
Retarded Time Equation Solutions 


Harmonic 

Number 

Case A 

OLD 

Case B 

NEW 

CaseC 

Case D 

Observer 1 ( x 1 

= .806 m, x 2 = 0.0 m 

, x 3 = -0.01m ) 


1 

142.63 

143.45 

142.59 

143.45 

2 

128.66 

130.61 

128.69 

130.61 

3 

117.54 

118.79 

117.62 

118.79 

4 

124.07 

124.86 

124.22 

124.86 

5 

116.24 

117.58 

116.16 

117.58 

6 

116.18 

114.90 

116.28 

114.90 

7 

106.50 

105.99 

108.44 

105.99 

8 

107.29 

102.56 

103.33 

102.56 


Observer 2 ( x 1 = 

.806 m, x 2 = 0.0m 

, x 3 = -.806 m ) 


1 

137.46 

138.29 

137.46 

138.29 

2 

129.65 

130.47 

129.65 

130.47 

3 

122.64 

123.44 

122.64 

123.44 

4 

116.38 

117.12 

116.38 

117.12 

5 

110.29 

110.99 

110.29 

110.99 

6 

104.99 

105.57 

104.99 

105.57 

7 

101.37 

101.83 

101.37 

101.83 

8 

96.72 

97.17 

96.72 

97.17 
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Figure 1.1 Example of the roots of the retarded time equation 


microphone 1 



<(> : wind tunnel microphone 
convention 

0 : ANOPP microphone 
convention 

r = 4 m 


Figure 2.1 ANOPP microphone convention and DNW microphone convention for the 
test case AN5 and N propeller. 
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SOUND PRESSURE LEVEL, dB 
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105 - 
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0 

Figure 3.1 
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MICROPHONE 1 
(0 = 90°, <> = 0°, r = 4m) 


M h = 0.872 
Cp = 0.049 


O DNW TUNNEL DATA 
4 NON-FLAT PLATE ANOPP PREDICTION 
A FLAT PLATE ANOPP PREDICTION 


T ■ 1 * ' 1 1” 

5 10 15 20 

HARMONIC NUMBER 


Noise comparison of ANOPP and DNW data for N propeller - test AN5. 
ANOPP predictions include the integral formulations for the boundary- 
layer thickness and the zero pressure gradient flat plate model. 
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SOUND PRESSURE LEVEL, dB 


120 f 


MICROPHONE 2 

(0 = 60°, <))=0 o > r = 4 m) 

Mu= 0.872 
Cp = 0.049 


O DNW TUNNEL DATA 
♦ NON-FLAT PLATE ANOPP PREDICT 
▲ FLAT PLATE ANOPP PREDICTION 


0 5 10 15 20 

HARMONIC NUMBER 


110 


• ft 


100 1 


t ° o 


90 A 


o o 


80 A 


* o ° 
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70 4 


60 4 


o 
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Figure 3.2 Noise comparison of ANOPP an DNW data for N propeller - test case AN5. 
ANOPP predictions include the integral formulations for the boundary- 
layer thickness and the zero pressure gradient flat plate model. 
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SOUND PRESSURE LEVEL, dB 


120 T 


115 - Cl 
x 
A 
* 

110 - ° 


105 A 


100 A 


95 A 


90 -f- 
0 


Figure 4.1 



HARMONIC NUMBER 


MICROPHONE 1 
(0 = 90°, $ = 0°, r = 4 m) 


□ DNW DATA 
A ICL=1 ,ICP=2 
A ICL=1,ICP=1 
O ICL=1,ICP=0 
+ ICL=0,ICP=2 
x ICL=0,ICP=1 
x ICL=0,ICP=0 


Noise comparison of different compressibility correction options of ANOPP 
predictions and DNW data for the test case AN5 and N propeller. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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SOUND PRESSURE LEVEL, dB 



MICROPHONE 2 

(0 = 60°, <J> = 0°, r = 4 m) 

□ DNW DATA 
A ICL=1,ICP=2 
A ICL=1,ICP=1 
O ICL=1 ,ICP=0 
+ ICL=0,ICP=2 
* ICL=0,ICP=1 
x ICL=0, ICP=0 


Figure 4.2 Noise comparison of different compressibility correction options of ANOPP 
predictions and DNW data from the test case AN5 for N propeller. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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OVERALL NOISE, dB 



Figure 4.3 Noise comparison for different compressibility options for the pressure coefficient. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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OVERALL NOISE, dB 



r = 4 m) 


ICL=0,ICP=2 

ICL=0,ICP=1 

ICL=0,ICP=0 


Figure 4.4 Noise comparison for different compressibility correction options for pressure 
coefficient - no correction applied for lift coefficient. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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OVERALL NOISE, dB 



HARMONIC NUMBER 


Figure 4.5 Noise comparison for different compressibility correction options for 
pressure coefficient. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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OVERALL NOISE, dB 



HARMONIC NUMBER 


Figure 4.6 Noise comparison for different compressibility correction options of 
pressure coefficient - Glauert theory applied for lift coefficient. 


ICL - Compressibility option for the lift coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 

ICP - Compressibility option for the pressure coefficient 
= 0 - no compressibility correction 
= 1 - Glauert compressibility correction 
= 2 - Karman-Tsien compressibility correction 
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Figure 5.1 Function f(x) for the case where the root is found in two iterations. 


F(X) = AX + BX - C + COS(X+D) 
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Figure 5.3 


Sound pressure level, dB 


Front view 
Diam. = 3.66 m 


Wind tunnel 

nozzle diameter = 7.32 m 



Figure 6.1 Microphone positions in 24-foot Wind Tunnel tests used for 
comparisons with the ANOPP predictions. 



Harmonic number 

Figure 6.2 Comparison of prediction with data for the 24-foot 

wind tunnel first test case (Dowty Rotol R212 propeller). 
( see Fig. 6.1 for microphone position ) 
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Sound pressure level, dB 



Harmonic number 

Figure 6.3 Comparison of prediction with data for the 24-foot wind 
tunnel second test case (Dowty Rotol R212 propeller). 

( see Fig 6.1 for microphone positions ) 



Sound pressure level, dB 



Harmonic number 

Figure 6.4 Comparison of prediction with data for the 24-foot wind 
tunnel third test case (Dowty Rotol R212 propeller). 

( see Fig. 6.1 for microphone position ) 
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ANOPP 60° 
position 


ANOPP 90° microphone 
microphone “t" position 



(j> = Wind Tunnel microphone convention 
0 = ANOPP prediction convention 


Figure 7.1 Microphone and nacelle positions in the DNW used for 
comparisons with the ANOPP predictions. 
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Sound pressure level, dB 


120 



Harmonic number 

Figure 7.2 Comparison of prediction and data from 90° directivity angle 
microphone for the DNW tunnel test. 
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Figure 7.3 Comparison of prediction and data from 90° directivity angle 
microphone for the DNW tunnel test. 
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Sound pressure level, dB 



Harmonic number 

Figure 8.1 Comparison of prediction with data from 60° directivity 
angle microphone for the DNW tunnel test. 
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Figure 9.1 Noise comparison of the first harmonic for grid combinations. 
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Figure 9.2 Noise comparison of the second harmonic for grid combinations. 
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igure 9.3 


Noise comparison of the third harmonic for grid combinations. 



Figure 10.1 Noise comparison of the first harmonic for grid combinations. 
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SOUND PRESSURE LEVEL, dB 




Figure 1 0.3 Noise comparison of third harmonic for grid combinations. 
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the propeller noise interior program. 
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LONGITUDINAL INDEX 
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PP code. 
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Figure 1 3.2 Amplitude of the first harmonic for the SR-7 ANOPP prediction. 
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Figure 13.4 Amplitude of the first harmonic of the SR-7 ANOPP prediction. 
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Figure 15.2 Amplitude of the first harmonic of ANOPP prediction for the UDF. 
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Figure 16.1 Source-to-observer positions used for the PTA propfan prediction 
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Figure 17.3 Predicted noise level of the PTA propfan. 
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